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Abstract 

Most of dynamic systems which exhibit chaotic behavior are also 
known to posses self-similarity and manifest strong fluctuations of all 
possible scales. The meaning of this terms is not always the same. In 
the present note we make an attempt to formulate the problem in the 
framework of functional analysis. The statistical hydrodynamics is 
taken vivid physical example. 

1 Introduction 

The problem of adequate mathematical description of hydrodynamic tur- 
bulence is one of the oldest but not yet solved problems in physics. The 
transition from laminar fluid flow to highly irregular chaotic regime was first 
discussed yet by Leonardo da Vinci. Later, when the Navier-Stokes equation 
(NSE) 

d t \ + v ■ Vv = -Vp + i/Av, V-v = 0, (1) 

which describes an incompressible fluid flow was written down, it was believed 
that ([!]) itself may contain all turbulence. It seems quite natural in the light 
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of modern developments in dynamical chaos, that even systems with a few 
degrees of freedom sometimes show unpredictable chaotic behavior. 

The turbulence problem is much more complicated. The velocity field 
v(x, t) has 3 x 3-D continuum degrees of freedom - thus we have a field 
theory problem. The dynamical object - the velocity field v(x, t) itself is a 
square-integrable function defined on IR 3 x IR 1 space. 

The crucial step in turbulence theory was done by A. N. Kolmogorov 
in a number of short papers |Kol41b| , |Kol41c| , [Kol41a|| , known as K41 the- 



ory. It was argued, that the turbulence, as a chaotic phenomenon, should be 
described in terms of random functions v(x, t, •) (The ideas of statistical de- 
scription of turbulence have been already discussed by Taylor [[Xay35fl ). It was 



suggested to consider the turbulence as a multi-scale phenomenon, formed of 
velocity fluctuations of all possible scales. The typical size L of dominat- 
ing fluctuations of average amplitude U is related to the Reynolds number 
Re = — , which determines the transition from laminar to turbulent flow. 
Following [[Fri95|1 we present K41 hypothesis in the form 



HI In the limit Re — > oo, all possible symmetries of the NSE, usually bro- 
ken by the mechanisms reducing the turbulent flow, are restored in a 
statistical sense at small scales and away from boundaries. 

H2 Under the same assumptions as in HI, the turbulent flow is self-similar 
at small scales, i.e.: 

5v{r, XI) '= X h 5v(r, I) (2) 

where v(r, XI) = v(r + XI) — v(r). (The equality-in-law means the 
coincidence of all statistical momenta.) 

H3 Under the same assumptions as in HI and H2, the turbulent flow has a 
finite non-vanishing rate e of dissipation of energy per unit of mass. 

Using these three assumptions Kolmogorov has derived the two-thirds law 

((5v(l)) 2 ) = Ce 2 ^ 3 , (3) 

the basic empirical law of fully developed turbulence. 

The scaling behavior (|3p of velocity field v means that we have to deal 
with functions of typical behavior 

\v(r + l) -v(r)\ ~l h (4) 
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with h < 1, ie. with non-differentiable functions; with nontrivial Holder 
exponent h = 1/3, if H3 hypothesis is strictly valid. The Holder condition 
\v(r + /) — v(r)\ < I is not very rare in physical problems: it can be found 
in condensed matter physics, quantum field theory etc. - - it becomes very 
significant only since the resolution (I) itself becomes a physical parameter, 
and that is why we ought to consider some objects of the form v(x, I, •), 
which are not yet defined. This problem is actively discussed amongst the 
turbulence community, but up to the authors knowledge is beyond functional 
analysis research, except for the problems directly related to wavelet analysis 

In present paper we investigate the problem of resolution dependent func- 
tions. The research was highly inspired by the fact, that there are at least 
two approaches to resolution dependent objects. First, based on the decom- 
position with respect to the group of affine transformations, is known as 
wavelet analysis. Second approach was exported from quantum field theory 
to turbulence and is based on the consideration of the Fourier transforms of 
functions cut at some large momentum k = K cut - ff- 

An attempt to consider the problem as a whole and join these two ap- 
proaches is presented below. 

2 Problem 

Self-similarity is a synonym of scale-invariance. To be scale-invariant means 
to have same properties at different scales. Classical fractals are scale- 
invariant by construction. Brownian motion is self-similar: if we look at 
the trajectory of Brownian particle at different resolution of a microscope we 
will observe more or less the same picture. 

As physical systems are considered, the word self-similarity is more fre- 
quently attributed to their dynamics than geometry. 

The self-similarity of hydro dynamic velocity field fluctuations {(Sv(l)) 2 ) ~ 
/ 2 / 3 is attributed to behavior of turbulent velocity field measured at different 
spatial scales. For hydrodynamical velocity field it is physically clear, that the 
measurement at scale Iq necessarily implies averaging of molecular velocities 
over certain space domain of typical size Iq. This procedure can be generalized 
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to "an averaging of a function up to scale 1" 



mm 



4n{x) = l- D [ 0(x-y)^y. (5) 

J\y\<l 

There are at least two conjectures here: 

1. The existence of "true" (with no scale) field 0;(x) : I — > 0. 

2. The homogeneity of the measure dfi(y) = d D y. 

Physically, it is quite clear, that two different fields </>z(x) and <?V(x) live in 
different functional spaces if I ^ V . It is meaningless, say, to subtract their 
values. Therefore, the velocity field of hydrodynamic turbulence is something 
more than a random vector field defined on R D x R. 

To characterize the turbulent velocity at certain point x we ought to know 
the collection of velocity values {0;(x)} at a set of scales labeled by I. The 
set of scales may be countable 



0: 



say k = 1/2 for period-doubling decomposition, or continuous. 

To characterize this set it was proposed in [ pG96|| to use a collection of 



unit fields at different scales - a "reference field" {TZi(x)}. The principal 
question arising here, is how to describe the interaction of the fluctuations of 
different scales. Practically, this problem is often coup with by decomposition 
of "real" (out-of- scale) field into slow (large-scale) and fast (small-scale) 
components. 

<p — V + v, where Mt> = 

(M means the averaging, or mathematical expectation, here and after.) In 
this approach the slow component V governs the equation for v and the 
even-order moments of v contribute to the equation for V. 

From the other hand, as we know from both Kolmogorov's theory and 
RG approach, there are no absolute scales in hydrodynamics, except for 
dissipative scale and external scale (the size of the system). So, at least 
at this middle - the Kolmogorov's range - the equations should be scale- 
covariant. The structure which reveals here looks like a fiber bundle over 
R D , with leaves labeled by scale. The fluctuations of different scales may be 
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dependent or independent for various physical situations; but at least some 
similarity should be present. 

To construct a basic system on this bundle let us follow the ideas of 
multi-resolution analysis [|Mal86 ] . Let us construct a system of functional 
subspaces {Vi : Vi C Ti.}, Where Ti is a space of physical observables. Let 
the system {Vi} be such, that 



1. 



c V 2 c Vi c V c 



2. 

nooVi = 0, U^v = n 

3. Subspaces Vj and Vj+i are similar in some way. 

f{x) e Vj <- f(kx) e Vs-u 
if {(pi(x)}i e i forms basis in Vj then {(f>i(kx)}i £ i forms basis in Vj-i 

If the sequence {V^} is bounded from above, the maximal subspace is called 
the highest resolution space. Let it be Vq. Then any function from V% can 
be represented as a linear span of V basic vectors. Therefore, the basis O 
of the highest resolution space provides a basis for a whole bundle. 

It seems attractive to generalize MRA axioms to the case of continuous 
set of scales. Since the chain of subspaces described above implies sequential 
coarse graining of the finest resolution field, some details are being lost in 
course of this process. The lost details can be stored into the set of orthogonal 
complements 

V = V 1 ® W u V l = V 2 @W 2} ... (6) 

So, Ti = J2k Wk, W}~r\Wj — if j k, and the system {Wi} can be considered 
instead of {V^}. The former has the structure of a- algebra, and thus it is 
suitable integration. 

The fact, that velocities (pi and <$>y live on different leaves suggest that 
their Fourier decomposition should be taken separately at their leaves 

fa(x) = J exp(-ik^x^)Uk (l Whk {l) ), (7) 
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or some other care should be taken about it in order not to mix fluctuations 
with the same wave vectors but contributing to different scales. The choice 
of the left-invariant measure dfi^(k^) is restricted by the fact, that velocity 
components measured at certain scale are mainly concentrated close to this 
scale. So the measure can be expressed as dfi^(k) = dkW(\l~ l — ak\), where 
W vanish at x — > ±00, a is a constant. 

The decomposition (0) turns to be a kind of Gabor transformation [ Gab46| . 
The measure can dfi can be explicitly scale-dependent, since the probability 
spaces (Qi,Ui,Pi) depend on scale.. 

At this point we arrive to the difference from standard wavelet approach, 
where the probability space is completely determined at finest resolution 
scale. However, if we accept the hypothesis that fluctuations of different 
sizes are statistically independent, we have to define the probability spaces 
separately. 



3 Wavelet realization 

We start our construction of multi-scale description with simplistic one- 
dimensional case, which is however of practical importance since only one 
component of velocity field is often measured. 

Any square integrable function f(t) G L 2 (R) can be represented as a 
decomposition with respect to the representations of affine group 

t r = at + b (8) 

o, s ^ 1 f 1 , ( t — b\ T _ r , dadb 

m = c ~< ■fvs' p {—) w * <a '" ) —' (9) 

which is just another form of the partition of unity with respect to represen- 
tation of affine group acting on a Hilbert space [ |GJP85|| 

1 = C^ f U*{g)\iP)d(, L (g){Tp\U{g), 

which holds if there exists such if) e Ti,, that 

CV := pn J G \^AU{g)^)\ 2 d^L{g) < 00 
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holds; dnL(g) denotes left-invariant measure on G. The scalar products 

W4g)f:= (f,U(9)*P) 

are known as wavelet coefficients. 

For the case of affine transformation group @ the normalization constant 
Gp can be easily evaluated in Fourier space 

c ^r\im ik=2 rm\ dK (10) 



where ip(t) — ^- j exp(zkt)ip(k)dk. For the affine group 



U(a,b)ijj{x) := -^=-0 ( - — - ) ; dfx L (a, b) = (11) 



The corresponding wavelet coefficients are 



W^a,b) = I (—)f(t)dt (12) 



a 



For practical analytical calculations it is often more efficient to calculate 
wavelet coefficients in Fourier representation, since multiplication should be 
done then instead of convolution: 

WJa,b)f = — ( y/aexp(ikb)ij(ak)f(k)dk (13) 

Z7T J 



and similar for reconstruction ([|). The decomposition (|T2| ) and its inverse 
(HD are known as wavelet analysis (See e.g. [Pau88|| for general review.) 
The scalar product ( jl2|) is readily seen to be the projection of the original 
"no-scale" function / to the subspace W a of MRA system |. 

If / is a random function defined on a probability space (fl,A, P), the 
wavelet coefficients 



W4a,b,-)= iJ-^t— )f(t,-)dt (14) 



a 



are also random; the stochastic integration is implied. As it is known from the 
theory of stochastic processes, any random function £(i,co>),t e IR 1 ,^ G Q 
can be represented in a spectral form 

£(*) = J^XMdX), (15) 
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where (pit, A) is a square-integrable function, T](dX) is a stochastic measure 

Mr](dX) = 0, M\7](d\)\ 2 = F(d\). 

The particular form of the spectral representation (|T5D is Fourier represen- 
tation 

f (t) = [ exp(i\t)rj(d\). 



In a multi-scale case we can introduce a collection of random processes, each 
of which belongs to its own leaf of MRA, labeled by resolution parameter a 

Ut) = J <t>a{t,\)r)a{d\). 

The peculiarity of stochastic case is, that in contrast to the decomposition 
of a function with respect to given basic wavelet ip(t), the function <f) a (t, A), 
which depends on both the properties of random process itself and filtering 
properties of measuring equipment, is not known exactly. Therefore, we have 
to construct a decomposition, which has a well defined limit to deterministic 
case and can be tackled without the exact specification of the form of basic 
wavelet ip. 

A straightforward way to do it, is to factorize the scaling part of the 
left-invariant measure from "purely stochastic part" 

£(t) = I Mt^)Va(dX) — 
J a 

or in the spectral form 

^(t) = 7^7T [ e^(Xa) Va (dX)-, (16) 

where rj a (d\) can be considered as generalized wavelet coefficients, the ex- 
istence of which does not require the existence of "no-scale" prototype. The 
left invariant measure dfi(a) = da /{2tt a) on the multiplicative group x' = ax 
instead of (|8]), since translations are already incorporated into the exponent; 
27r multiplier is introduced for the convenience of Fourier transform. 



The representation (|lj]) was constructed only to meet the non-stochastic 
limit and is not unique. For instance, we can redefine the spectral measure to 
incorporate both the properties of the signal and that of measuring apparatus 

m = I e^%(dX)-. 
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The specific energy-per-scale density can be easily evaluated then 



JE(a)da = Mf£{t)Z(t)dt 

= M/exp(zt(Ai - \ 2 ))r] ax {d\x)ri a2 {d\2)d^{a 1 )dii{a2)dt , , 

= M J r] ai (d\)fj a2 (d\)dfi(a) 

= J F(ai, a<2\ d\)dfi(ai)dfi(a 2 ) = E 

The equation QT7| ) is a stochastic counterpart of a well known equation for 
wavelet energy per scale 

^Hc-ZlM* i /*f it (18) 

Similar spectral characteristics have been already used for the analysis of 
turbulent data ||Ast96 . 



If random functions are considered the equation (TPS!) can be rewritten in 
RG-like form 

dE 1 



d In a 2tiC., 



M J \W(a,k)\ 2 dk, (19) 



where W(a, k) = ip(ak)f(k) can be understood as the original noisy signal 
/ perceived by filter ip, i.e. as generalized wavelet coefficients, not necessary 



having "no-scale" prototypes, cf.eq.|16 



The logarithmic derivative at the 1. h. s. of fllTf ) is exactly that of renor- 
malization group equation. However, the r. h. s of this equation was obtained 
without any cutoff assumptions, it was formally derived from the decompo- 
sition of initial ("infinite resolution") signal filtered by window function ip. 

Physically, the existence of the infinite resolution limit is often meaning- 
less. For instance, the behavior of electromagnetic coupling constant at Plank 
scale is just a nonsense. The same happens in hydrodynamics when the scale 
parameter comes close to the mean free path. However, there is a principal 
difference between hydrodynamical turbulence description and quantum field 
theory. The scale / (resolution) becomes a physical measurable parameter in 
hydrodynamics, and thus W(a, k) can be considered as Fourier components 
of the velocity field fluctuations of different typical sizes. The study of the 
behavior of W(a, k) can provide more consistent picture of what happens at 
different scales than standard Fourier decomposition. 



9 



4 High frequency cut-off 



Being utmost scale-invariant at moderate scales, the behavior of turbulent 
velocity field changes when approaching the smallest and largest scales, be- 
tween which the hydrodynamical description is valid. The former is the 
Kolmogorov dissipative scale (rj), the latter is the size of the system. The 
size of the system can often be set to infinity with no harm to physics; whilst 
the dissipative scale is of physical importance, since the energy dissipation 
rate e is very constant which determines the turbulence behavior in inertial 
range. 

That is why in RG, as well as in spectral calculations, the cut-off depen- 
dent velocity field is often considered 



The cumulative energy of all harmonics with wave vectors less or equal to 
the cut-off value F is one of the main spectral characteristics of developed 
turbulence 



Similarly, we can consider the cumulative energy of all velocity fluctuations 
with typical size greater or equal to a given A. For simplicity let us consider 
a one-component velocity field considered as a function of time 




(20) 




(21) 



E(A) 




2 dadb 



a- 



(22) 



where 




and 




is the total energy of all velocity fluctuations. 
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For non- vanishing A 



roo \ip(y)\ 2 i 

E(A) = F(A)E, where, F(A) = — - . (23) 

For a better definiteness, let us calculate the filtering function F(A) for 
a particular family of vanishing momenta wavelets 

^ n {x) = (-l) n — exp(-x 2 /2), Mk) = V2^(-ik) n exp(-k 2 /2), (24) 
ax n 



often used for studying of hydrodynamical velocity field ||MBA91| , [MBA93 
The normalization constant for this family is 

/oo 
k 2n ~ x e- k dx = 2nT(n) 
-oo 

and so 

K{A) = /sjgp. (25) 

The derivative of cumulative energy with respect to logarithmic measure 
da /a is 

9E _ E 9Jm = M # 



dlnA d\nA d\nA 
where f n (x) = x n ~ 1 e~ x /T{n). So we arrive at RG like equation 

dE 2A 2n exp(-A 2 ^ 



dlnA T{n) 



E (26) 



For sufficiently small A the exponential term is close to unity, and thus the 
behavior is approximately proportional to A 2 . 



5 Self-similarity in a bounded domain of scales 
and Scale-relativity 

A direct approach to scale dependent functions, which was proposed by 
Nottale [fSlot93| , is based on the assumption that physics is scale dependent, 
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but scale covariant. The former means to give up the differentiability and 

consider functions v (x, I) with lim — — - — not necessary existing; the latter 

i->o ox 

imposes generalized scale covariant equation 

^T = ^M (27) 

which states that scale behavior of a scale-dependent function is completely 
determined by the value of this function at a given scale. 

The benefit of scale relativity |Not9^] approach is to provide a possibility 



to account for the processes which admit self-similarity only at a limited 
domain of scales 

v(Xl) = \ h v(l),r] <I<A. 

For a homogeneous scale dependent function (mono-fractal) v = vo(X/r) 5 
where 5 is independent of scale r. To get the additive, rather than multi- 
plicative form for two sequential scale transformations, it is convenient to 
express them in logarithmic form : 

where 5(r') = 8(r) = 5, V = ln(r/V). In real hydrodynamic turbulence the 
multi-fractal behavior is observed 

5 = 6{r). 

In the inertial range ?| C r < A the exponent 5 is practically a constant, 
but closer to the limiting scales rj and A the dependence of scale becomes 
significant. 

The logarithmic form suggests a direct generalization on the multi- 
fractal case 5 = 5(r), which is like the generalization of Galilelian transform 
to Lorenzian transform: 

X' = T(¥)[X-¥T] 

V = T(V)[A(V)X + B(V)T], { > 



where T = \n(l/l ),X(T) = lnM(uj/u ), see jDG9| for details. 



12 



The composition of two scale transformations of this type behaves like 
a composition of two Lorentzian boosts: when approaching the unpassible 
limit, then instead of "Galilelian" 

5(r) =S (1- ln 2 (A/r)/ ln 2 (A/^))- 1/2 (30) 



for the power behavior of velocity field. The transformation group fl29|) can 
be used to construct a wavelet decomposition at a limited domain of scales. 
This, however, is the subject of the next paper (in preparation). 



6 Conclusion 

In present paper we give a mathematical framework for the analysis of func- 
tions which depend on scale. Usually, the scale-dependent functions express 
the value of a certain physical quantity measured at a point x by averaging 
over a box of size I centered at x. Such functions are often used in hydrody- 
namics, geophysics, signal analysis. One of the most known ways of treating 
resolution-dependent functions is to identify the size of the box (/) with the 
inverse wave number of Fourier transform k" 1 = I. The higher wave num- 
bers are then cutted off. Sometimes this procedure leads to confusion (many 
problems of field theory approach to hydrodynamic turbulence originate from 
this confusion). 

In our approach, using the ideas of wavelet analysis, we keep wave vectors 
(k) and scales (a) separately. We derive a renormalization group like equa- 
tions, which can be used to study the energy distribution between different 
scales. The results can be applied for the investigation of turbulence velocity 

experimental data, as well for further theoretical research. 
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